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SECTION  I 


INTRODUCTION 

The  aerodynamicists  ability  to  analyze  transonic  flows  has  advanced 
significantly  in  the  past  decade.  Transonic  airfoil  codes  are  now  as 
accurate  as  subsonic  airfoil  codes  and  transonic  wing  codes  are  rapidly 
approaching  the  accuracy  of  subsonic  wing  codes.  Currently  effort  is 
being  directed  toward  extending  wing  codes  to  include  body  effects  and  the 
effects  of  other  components. 

These  transonic  codes  are  much  more  costly  to  run  than  their  corres¬ 
ponding  subsonic  codes  because  of  their  non-linear  formulation.  This  extra 
cost  requires  that  specific  benefits  of  using  these  codes  must  be  identi¬ 
fied.  One  way  to  ascertain  the  value  of  these  advances  is  to  apply  these 
codes  to  the  redesign  of  wings  that  were  designed  prior  to  the  development 
of  these  transonic  codes.  The  goal  of  such  redesign  efforts  would  be  to 
quantify  the  improvement  in  the  transonic  cruise  efficiency  of  a  given 
aircraft  due  to  the  use  of  the  transonic  codes. 

A  prime  candidate  for  such  a  redesign  effort  is  the  C-141  aircraft. 

A  large  part  of  its  flying  time  is  spent  in  cruising  at  transonic  speeds, 
approximately  M=0.74,  and  since  the  C-141  fleet  consumes  about  15  percent 
of  the  total  Air  Force’s  jet  fuel  annually,  a  small  improvement  in  cruise 
efficiency  could  net  a  significant  savings  in  fuel. 

The  Flight  Dynamics  Laboratory’s,  FDL,  interest  in  transonic  codes 
and  their  application  to  the  redesign  of  an  existing  aircraft’s  wing 
merged  nicely  with  Lockheed -Georgia ’ s  ongoing  effort  to  reduce  the  creep 
drag  of  the  C-141  aircraft.  Under  IRAD  work  they  had  isolated  the  creep 
drag  problem  of  the  C-141  to  the  airfoil  profile.  They  had  also  done 
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some  basic  cut  and  try  experimental  work  to  determine  the  effect  of  the 
leading  edge  shape  on  the  airfoil’s  creep  drag.  Based  on  this  work  FDL 
and  Lockheed-Georgia  formulated  an  effort  to  use  transonic  airfoil  and 
wing  codes  to  optimize  the  0141  wing. 

Under  contract  to  FDL,  Lockheed-Georgia  produced  leading  edge  modifica¬ 
tions  to  the  0141  wing  to  reduce  its  cruise  drag.  These  modifications  were 
restricted  to  the  first  twelve  percent  chord  due  to  the  position  of  the 
wing’s  front  spar.  The  two  resulting  designs  were  tested  at  Arnold  Engineer¬ 
ing  and  Development  Center  (AEDC)  16-foot  Transonic  Tunnel  and  the  results 
are  documented  in  Reference  5.  Reference  5  also  contains  background  on 
these  efforts,  a  description  of  the  procedures  used  to  design  the  wing  leading 
edge  modifications,  and  a  description  of  two  other  wing  modifications  tested, 
a  wing  swept  tip  extension  and  trailing  edge  anti-drag  bodies.  Due  to  in¬ 
accuracies  in  model  construction,  the  test  results  were  not  totally  satis¬ 
factory  and  FDL  and  Lockheed-Georgia  entered  into  a  cooperative  program  to 
design,  build,  and  test  a  third  leading  edge  modification.  That  effort  and 
its  results  are  the  subject  of  this  technical  report. 
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SECTION  II 


AERODYNAMIC  DESIGN 

1.  BACKGROUND 

The  C-141  exhibits  an  increase  in  drag  with  increasing  Mach  number 
throughout  the  subsonic  Mach  number  range  in  wind  tunnel  tests  and  in 
flight.  Figure  1  shows  this  creep  drag  and  drag  rise  for  the  C-141A. 
Figure  1  also  shows  that  a  subsonic  drag  rise,  of  the  same  order  of  magni¬ 
tude  as  that  of  the  airplane,  has  been  measured  on  the  wing  alone  and  on 
the  basic  airfoil. 

a.  TWO-DIMENSIONAL  DESIGN  AND  TEST 

Previous  studies  by  Lockheed-Georgia  have  attempted  to  isolate 
the  source  of  the  C-141  creep  drag  and  drag  rise.  These  studies  indicated 
that  the  problem  was  not  related  to  either  induced  drag  changes  or  pre¬ 
mature  flow  separation  as  Mach  number  increases.  Based  on  these  results 
they  then  focused  on  the  more  likely  source  of  increased  profile  drag  due 
to  viscous  and  pressure  losses.  IRAD  studies  were  initiated  at  Lockheed- 
Georgia  in  an  effort  to  understand  more  fully  the  relationships  between 
pressure  distributions  and  the  drag  rise  phenomena  as  they  occur  on  2-D 
airfoils  such  as  that  of  the  C-141A. 

Initial  efforts  compared  the  drag  rise  characteristics  of  con¬ 
ventional  airfoil  shapes  such  as  the  C-141  airfoil  to  a  number  of  super¬ 
critical  airfoils.  Two  geometric  characteristics  emerged  as  important 
factors  in  airfoil  drag  rise  performance.  They  were  the  trailing  edge 
cusp  region  and  leading  edge  contour.  An  example  of  the  effect  of  leading 
edge  contour  is  shown  in  Figure  2.  The  two  research  airfoils  in  this  test 
were  essentially  the  same  except  for  the  forward  15  percent  of  the  chord. 
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Moving  the  leading  edge  of  the  airfoil  upward  results  in  a  moderate  pres¬ 
sure  peak  on  the  first  10  percent  of  the  upper  surface  with  a  weaker  shock 
development  which  is  further  aft  than  the  shock  on  the  original  airfoil. 
This  influence  of  the  leading  edge  suction  (or  pressure)  peak  is  evident 
throughout  the  Mach  number  range,  causing  improved  creep  drag  as  well  as 
delayed  drag  divergence  performance. 

Based  on  these  results,  Lockheed  expanded  its  IRAD  studies  to 
include  the  design,  test,  and  analysis  of  a  series  of  leading  edge  modifi¬ 
cations  to  a  two-dimensional  C-141  airfoil.  For  the  C-141,  the  location 
of  the  wing  spar  restricted  the  extent  of  modification  to  the  forward  12% 
chord,  thus  limiting  the  amount  of  drag  reduction  which  might  be  achieved. 
An  existing  2-D  model  of  a  C-141  airfoil  at  r)  =  0.389  was  cut  at  the  12% 
chord  location  and  modified  to  accept  removable  leading  edges.  To  provide 
a  baseline  for  comparison  a  replacement  leading  edge  was  constructed  using 
the  original  airfoil  ordinates. 

Manual  iteration  of  the  viscous  airfoil  theory  of  Bauer,  et  al^^ 

was  used  to  develop  a  number  of  designs.  Also,  one  design  was  determined 

during  a  cooperative  effort  between  Lockheed-Georgia  and  NASA/Ames  using 

(2) (3) 

the  2-D  CONMIN  method  of  optimization.  This  latter  design  procedure 

was  also  used  to  verify  the  previous  manual  designs. 

Results  for  the  baseline  and  two  of  the  modified  leading  edges 
are  summarized  in  Figure  3.  The  2-D  CONMIN  leading  edge,  LE6,  uses  the 
existing  lower  surface  shape  and  is  thus  a  more  practical  modification* 
Several  larger  per tubations ,  such  as  LE  3A,  were  included  to  provide  a 
range  of  leading  edge  shapes.  Figure  3(b)  illustrates  relative  changes  of 
the  flow  over  the  leading  edges  and  demonstrates  that  weaker  shock  forma¬ 
tions  were  achieved  with  the  modifications.  The  drag  rise  data  of 


Figure  3(c)  shows  the  substantial  improvement  of  both  LE6  and  3A  over  the 
baseline.  Although  the  expected  reductions  in  creep  drag  at  the  low  and 
intermediate  Mach  numbers  (M<0.65)  did  not  materialize,  creep  drag  is 
reduced  for  M>0,65  and  an  increase  in  drag  divergence  Mach  number  of 
approximately  0.02  is  achieved.  At  M2D=0.73,  corresponding  to  a  cruise 
Mach  number  of  0.775,  a  drag  reduction  of  11  counts  was  indicated.  LE6 
proved  to  be  as  good  or  beLter  than  any  of  the  others,  thus  verifying  the 
capability  of  the  CONMIN  design  approach  and  also  emphasizing  that  the 
modification  may  be  applied  to  the  upper  surface  only. 

At  the  section  lift  coefficients,  which  correspond  to  the  equiva¬ 
lent  2-D  local  lift  coefficients  at  cruise  on  the  mid  and  outer  span  por¬ 
tions  of  the  C-141  wing,  drag  reductions  of  11  to  16  counts  were  typical. 
Since  15  counts  of  drag  represents  about  6  percent  of  the  total  C-141  cruise 
drag,  a  design  goal  of  5  to  7  percent  drag  reduction  at  cruise  for  the  wing 
was  adopted  as  a  standard  for  the  3-D  wing  aerodynamic  design, 
b.  PRIOR  THREE-DIMENSIONAL  DESIGN  AND  TEST 

The  two-dimensional  results  discussed  above  were  presented  to 
Flight  Dynamics  Laboratory  personnel.  This  application  of  transonic  aero¬ 
dynamic  codes  to  improving  on  existing  airfoil  and  wing  designs  fit  nicely 
into  FIMMs  External  Aerodynamics  Group’s  charter.  Thus  an  effort  was 
contracted  for  Lockheed-Georgia  to  extend  their  2-D  work  to  a  3-D  design 
that  would  be  tested  at  AEDC’s  16T  wind  tunnel. 

Three  leading  edge  modifications  were  designed  from  which  two  were 
selected  for  model  fabrication  and  testing  in  conjunction  with  the  basic 
C-141  leading  edge.  The  first  modification  was  designed  using  a  Lockheed- 
Georgia  version  of  2-D  CONMIN  extended  to  permit  a  wider  range  of  design 
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variables.  A  3-D  numerical  optimization  scheme  under  development  by  Hicks 
similarly  extended  in  scope,  was  employed  to  obtain  a  second  leading  edge. 

The  third  was  based  on  the  best  of  the  2-D  modifications  tested  in  the 
Lockheed-Georgia  CFWT.  The  2-D  and  3-D  CONMIN  designs  showed  close  similari¬ 
ties  in  ordinates,  theoretical  pressure  distributions  and  drag  coefficients; 
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therefore,  only  the  3-D  CONMIN  wing,  W  ,  was  tested  and  not  the  2-D.  The 

36 

second  wing  tested,  W  ,  was  based  on  the  earlier  2-D  work. 

An  existing  model  of  the  C-141B,  0.044  scale,  was  modified  for  this 
test.  The  wing  was  cut  at  12%  chord  to  remove  the  leading  edge.  Three  new 

sets  of  leading  edges  were  fabricated,  one  set  to  replace  the  baseline  C-141 
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leading  edge,  W  ,  and  two  modified  sets,  W  and  W  ,  for  evaluation. 

The  results  of  this  test,  as  documented  in  References  5,  6,  and  7, 

were  disappointing  in  terms  of  drag  reduction  measured.  The  best  leading 
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edge,  according  to  the  AEDC  test  data,  was  the  W  for  which  a  drag  reduction 
at  cruise  of  approximately  5  to  9  counts  (2.0  to  3.6%)  was  obtained,  Figure  4 
Measurement  of  the  model  after  the  test  revealed  inaccuracies  in  the  con¬ 
struction  of  the  basic  and  modified  leading  edges,  Figure  5.  These  errors 

35 

were  such  that  the  geometric  differences  between  the  manufactured  W  and 
12 

W  were  much  less  than  they  should  have  been.  Subsequent  analysis  of  these 
measured  ordinates  by  the  transonic  codes  used  in  the  design  process  revealed 
drag  levels  in  excellent  agreement  with  the  test  results.  Based  on  this 
correlation  it  was  concluded  that  the  target  drag  reduction  goal  of  between 
5  and  7  percent  would  have  been  obtained  with  the  proper  leading  edge 
ordinates. 

This  effort  was  successful  In  demonstrating  the  capabilities  of 
the  transonic  codes  to  modify  an  existing  wing  to  obtain  a  cruise  drag 
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reduction.  The  target  drag  reduction  of  6%  was  not  achieved  due  to  model 
manufacturing  inaccuracies.  However,  the  increase  in  leading  edge  suction 
did  occur  and  this  net  thrusting  effect  tended  to  reduce  creep  drag  and 
increase  the  drag  divergence  Mach  number  as  predicted. 

The  results  were  promising  enough  that  the  Flight  Dynamics  Labora¬ 
tory  and  Lockheed-Georgia  decided  to  continue  the  effort  on  a  cooperative 
basis.  A  new  leading  edge  would  be  designed  and  built  based  upon  the 
knowledge  gained  from  this  effort.  FDL  would  then  provide  the  funding  for 
an  AEDC  16T  wind  tunnel  test.  This  effort  is  documented  in  this  report, 

2.  DESIGN  METHODOLOGY 

Four  different  design  methodologies  have  been  used  during  this  effort. 
These  methodologies  were  used  to  evaluate  the  relative  benefits  and  penal¬ 
ties  of  using  each.  Three  (2-D  CONMIN,  2-D  Airfoil,  and  3-D  CONMIN)  were 
used  in  designing  the  leading  edges  for  the  first  three-dimensional  test 
at  AEDC.  The  fourth  method  (2-D  equivalent)  was  used  to  design  the  leading 
edges  for  this  entry  at  AEDC. 

a.  PRIOR  DESIGN  APPROACHES 

For  the  first  entry  three  leading  edges  were  designed  but  only 
two  tested.  The  leading  edge  designed  by  using  the  2-D  CONMIN  method  was 
so  close  geometrically  to  the  3-D  CONMIN  leading  edge  that  it  was  not  built 
or  tested.  The  other  leading  edge  tested  was  based  on  the  best  of  the 
earlier  2-D  airfoils  from  the  Lockheed  IRAD.  A  summary  of  these  design 
procedures  is  shown  in  block  diagram  form  in  Figure  6.  Details  of  the 
individual  design  methods  are  given  in  References  6  and  7. 

b.  2-D  EQUIVALENT  APPROACH 

This  approach  was  developed  after  the  AEDC  16T  wind  tunnel  test  in 
April  1978.  It  was  based  on  the  experience  gained  in  using  the  three  design 
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approaches  mentioned  earlier.  Several  important  lessons  were  learned  or 
reinforced  from  these  previous  design  and  test  efforts.  These  lessons 
include: 

(1)  2-D  and  3-D  numerical  optimization  give  similar  results  for 
high-aspect-ratio,  moderately  swept  clean  wings. 

(2)  Numerical  optimization  using  inviscid  theory  produced  a  de¬ 
sign  which  performed  as  well  as  a  more  costly  design  based  on  viscous 
theory.  This  result  can  be  expected  where  the  design  changes  do  not 
significantly  affect  the  boundary  layer  properties. 

(3)  Very  tight  control  of  r-o del  fabrication  tolerances  is  re¬ 
quired  when  minor  leading  edge  cont c<;  .hanges  are  being  evaluated. 

The  first  two  lessons  above  combined  to  make  a  2-D  equivalent, 
inviscid  approach  in  the  optimization  process  desirable  in  terms  of  holding 
down  computer  costs.  Use  of  a  3-D  viscous  code  was  limited  to  analysis  of 
the  final  wing  design.  The  third  lesson  impacted  the  model  manufacturing 
process  which  is  discussed  in  Section  III. 2. 

The  equivalent  2-D  airfoil  at  each  control  station  was  obtained 
using  simple  sweep  theory  to  convert  from  3-D  to  2-D  and  vice  versa.  The 
equivalent  2-D  airfoils  and  conditions  are  defined  by: 

(y/c)^D  ^Y^c^treamwise  ^  cos 

M2D  ^streamwise  X  C0S  ^EFF 

C  =  C  X  - r-i - 

P2T)  P3D  cos2 

where  Agpp  ■  Sweep  at  Upper  Surface  Shock 
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For  this  design  effort  the  number  of  control  stations  was  five;  r|  =  0.118, 
0.445,  0.625,  0.815,  and  0.95.  The  root  station,  o  =  0.118,  was  held  as 
the  unmodified  basic  airfoil  and  the  other  four  were  optimized. 

The  2-D  CONMIN^  code  was  again  used  for  the  optimization  but  it 

(9) 

was  linked  to  Jameson’s  full  potential  inviscid  transonic  airfoil  code 
rather  than  the  Bauer  code  used  in  the  earlier  2-D  C0NM1N  approach.  The 
function  to  be  minimized  was  wave  drag  at  each  span  station.  The  final 
design  produced  by  this  2-D  equivalent  approach  was  analyzed  using  a 
viscous  version  of  the  Bauer  code^^ . 

The  choice  of  design  conditions  was  a  critical  step  for  this  effort, 
as  the  leading  edge  shape  was  found  to  be  very  sensitive  to  design  lift 
coefficient.  The  design  cruise  Mach  number  for  this  effort  was  selected 
as  M  -  0.79  for  two  reasons.  First,  this  Mach  number  provides  a  strong 
shock  wave  and  hence  greater  resolution  in  the  sensitive  calculations  made 
during  the  minimization  process.  Second,  the  cruise  range  parameter, 

M(L/D) ,  from  the  previous  AEDC  test  of  the  modified  leading  edge  was  optimum 
at  M  5  0.79.  The  appropriate  2-D  lift  coefficient  for  the  transonic  program 
was  selected  by  matching  the  upper  surface  leading-edge  pressures  with  the 
previous  AEDC  test  data  at  each  control  station.  Because  of  the  higher  than 
expected  loss  in  lower  surface  loading  due  to  fuselage/wheel-pod  overpres¬ 
sures  measured  in  the  AEDC  test,  the  2-D  section  lift  coefficients  are 
somewhat  higher  than  those  used  in  the  previous  2-D  design. 

The  leading  edge  shape  was  defined  as  a  four th-order  polynomial, 

Bf  B2  B^ 

y  »  A  +  A,  Y  +A„y  +  A~  Y  +  A,  y  where  the  A’s  and  Bfs  are  in- 
o  1  z  J  4 

dependent  design  variables.  Both  ordinate  and  slope  continuity  were  en¬ 
forced  at  the  match  points  on  the  wing  box.  A^  and  A^  were  selected 
to  satisfy  these  match  conditions,  thus  leaving  seven  design  variables. 
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Additionally,  angle  of  attack,  a,  was  used  as  a  "dummy”  design  variable 
to  provide  the  required  lift  coefficient.  Starting  values  for  the  Afs  and 
Bfs  were  determined  by  curve  fitting  the  basic  C-141  leading  edge.  The 
airfoil  chord  was  held  fixed. 

A  great  amount  of  detail  design  and  analysis  was  done  with  the 

r)  «  0.794  station  and  the  results,  summarized  in  Figure  7,  demonstrate 

the  important  aspects  of  this  design  effort.  Figure  7  is  a  convenient 

means  of  visualizing  the  amount  of  change  in  ordinate  at  2-percent  chord 

and  its  corresponding  effect  on  wave  drag.  Two  types  of  theoretical  wave 

drag  are  shown,  one  for  a  viscous  code  analysis  and  the  other  for  an  in- 

viscid  code.  The  similarity  in  the  curves  validates  the  use  of  an  inviscid 

code  in  the  optimizer.  The  first  three  data  points  shown  are  for  airfoils 
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from  the  basic  wing,  W  ,  the  3-D  CONMIN  wing,  W  ,  and  a  wing  designed 

39  40  41 

using  this  2-D  equivalent  approach,  W  .  Airfoils,  W  and  W  ,  with 

39 

greater  amounts  of  modification  than  the  optimized  airfoil,  W  ,  were 
determined  by  a  manual  iteration  of  the  curve  fitting  process.  These  re¬ 
sults  substantiate  the  existance  of  a  drag  bucket  and  show  the  eventual 
drag  rise  as  a  reflex  in  surface  curvature  on  the  upper  surface  occurs  due 

to  the  12  percent  chord  match  point.  Relaxation  of  this  match  point  con- 

37 

straint  to  17  percent  chord  produced  the  W  airfoil,  which  is  a  substan¬ 
tial  improvement  over  the  earlier  designs. 

Relaxation  of  the  match  point  to  17  percent  chord  was  also  used 
at  the  r|  *  0.95  station  after  initial  modification  efforts  using  the  12% 
match  point  failed  to  produce  a  drag  reduction.  At  r\  =  .625  and  n  ”  *445 
the  match  point  was  moved  back  to  12  percent  chord.  At  the  r|  =  .445 
station,  the  optimized  airfoil  differed  little  from  the  basic  airfoil  and 
offered  no  aerodynamic  improvement  over  it.  The  upper  surface  contour  at 
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n  *  0,794  was  duplicated  at  the  n  *  0.445  station  and  the  analysis  program 

37 

showed  the  Improvements  noted  in  Figure  8  for  the  W  configuration. 

Figures  9  and  10  show  the  results  at  n  *  0.625  and  r)  *  0.95.  The  basic, 
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W  ,  airfoils  and  modified,  W  ,  airfoils  are  compared  in  Figure  11.  These 
airfoils  were  selected  for  testing  in  this  effort. 

The  theoretical  pressure  distributions  for  the  basic  and  modified 
airfoils  are  shown  in  Figure  12,  These  predictions  are  from  the  Bauer ^ 
program.  Comparisons  of  these  pressure  distributions  reveal  remarkable 
improvements  in  shock  formation  considering  the  limited  extent  of  the 
modification.  There  are  also  significant  changes  in  the  upper  surface  pres¬ 
sures  before  the  shock.  The  major  change  is  the  creation  of  a  pressure  peak 
on  the  first  12%  of  the  upper  surface  but  the  influence  of  this  peak  does 
carry  over  all  the  way  to  the  shock.  On  the  lower  surface  the  change  in 
pressure  coefficients  is  much  less  with  the  largest  effect  being  at  the  two 
outboard  stations, 
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A  comparison  of  W  to  W  shows  the  difference  in  results  between 

the  two  optimizations  and  the  basic  wing.  Figure  13  summarizes  the  span- 

wise  variation  of  the  y/c  ordinate  at  2  percent  chord  for  the  three  wings. 

The  increasing  deflection  spanwise  for  the  modified  wings  effectively 

removes  some  of  the  existing  large  leading-edge  camber  as  the  tip  is 

approached.  Figure  14  shows  the  resulting  wave  drag  reductions  across 

the  wing  span  for  these  modif ications.  Integration  of  these  curves  shows 
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a  reduction  in  wave  drag  of  8  counts  for  W  and  5  counts  for  the  VT  de- 

35 

sign.  In  the  AEDC  test  of  the  incorrectly  fabricated  W  the  measured  drag 
reduction  was  about  9  counts  which  indicates  that  the  3-D  drag  improve¬ 
ments  may  be  substantially  higher  than  indicated  in  Figure  14.  Thus  the 


11 


37 

expected  drag  reduction  for  W  was  approximated  at  14  to  16  counts  at 

M  =  0.79.  The  primary  reasons  for  the  greater  leading  edge  deflection  and 
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a  corresponding  greater  reduction  in  drag  for  W  versus  W  are  relaxation 
of  the  12  percent  constraint  on  the  outboard  wing  and  use  of  higher  2-D 
design  lift  coefficients. 

3.  ADDITIONAL  DRAG  REDUCTION  CONCEPTS 

Two  additional  wing  modifications  from  Lockheed-Georgia f s  IRAD  were 
included  in  the  wind  tunnel  test  portion  of  this  study  to  evaluate  addition¬ 
al  drag  reduction  potential  for  the  C-141  aircraft.  These  two  modifications 
were:  swept  wing  tips  (extended  span),  and  a  series  of  wing  trailing  edge 
bodies  called  anti-drag  bodies  spaced  across  the  wing.  The  necessary  model 
components  for  these  configurations  were  provided  at  the  expense  of  the 
Lockheed-Georgia  Company  and  the  test  time  for  the  earlier  AEDC  entry  was 
included  in  the  FDL  funded  test  program.  For  this  entry  Lockheed-Georgia 
Company  funded  the  test  time  for  these  items  while  FDL  funded  the  leading 
edge  test  time. 

During  the  earlier  test  entry  these  two  concepts  were  tested  individual¬ 
ly  on  the  basic  wing.  For  this  test  they  were  tested  together  on  the  modi¬ 
fied  wing  to  see  if  the  drag  savings  from  the  individual  concepts  were 
additive. 

a.  SWEPT  WING  TIPS 

There  are  two  reasons  for  trying  swept  wing  tips  on  the  C-141. 

First,  the  obvious  increase  in  aspect  ratio  reduces  the  induced  drag  for 
a  given  wing  lift  coefficient  while  the  wing  loading  reduction  for  a  given 
total  lift  should  reduce  shock  losses  on  the  wing.  Second,  sweeping  of  the 
tip  should  counter  the  tendency  for  the  isobars  on  the  present  wing  to 


become  unswept  at  the  tip.  The  amount  of  extension  and  sweep  for  this 
modification  was  based  on  experience  with  a  similar  application  on  the  C-5A 
wing.  The  planform  of  the  swept  wing  tip  extension  is  shown  in  Figure  15. 
b.  WING  ANTI-DRAG  BODIES 

Lockheed-Georgia  experience  with  development  tests  on  the  C-5A 
flap-track  fairings  indicated  that  such  fairings  could  be  designed  to  pro¬ 
vide  a  net  drag  reduction  notwithstanding  the  increased  profile  drag  due  to 
the  fairings  themselves.  Analysis  of  the  C-5A  results  show  that  the  major 
effect  at  cruise  speeds  must  be  a  reduction  in  shock  strength  due  to  an 
effective  change  in  camber.  For  application  to  the  C-141  wing,  an  approach 
to  design  a  set  of  anti-drag  bodies,  similar  to  flap-track  fairings,  was 
devised . 

Assuming  that  a  i ^vorable  change  in  2-D  camber  could  be  applied  to 
the  3-D  wing  by  means  of  properly  sized  anti-drag  bodies,  use  was  made  of 
the  2-D  Bauer^  airfoil  analysis  program  to  determine  the  effects  of 
changes  in  camber.  To  increase  the  aft  camber  in  a  manner  which  could  be 
approximately  represented  by  isolated  bodies  at  the  trailing  edge,  a  lower 
surface  modification  was  devised  as  shown  in  Figure  16(a) .  The  lower  sur¬ 
face  ordinates  were  modified  by  the  addition  of  thickness  which  varied  from 
zero  at  X/C  =  0.65  to  2  percent  at  the  trailing  edge.  Figure  16(b)  shows 
the  effect  on  the  theoretical  pressure  distribution.  In  comparison  with 
the  basic  airfoil,  the  higher  aft  loading  permits  a  reduction  in  section 
angle  of  attack  from  2.7  degrees  to  0.9  degrees.  This  results  in  a  de¬ 
crease  in  upper  surface  shock  strength  coupled  with  a  rearward  shock  move¬ 
ment,  which  reduces  the  section  compressibility  drag  from  24  counts  to 
5  counts. 
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The  increase  in  lower  surface  trailing  edge  thickness  of  2  percent 


chord  was  converted  into  a  spanwise  cross-sectional  area  distribution  at 
the  trailing  edge.  To  ensure  a  reasonable  representation  of  this  distribu¬ 
tion  with  adequate  spacing  to  prevent  local  channel  flow  interference 
effects,  a  total  of  eight  bodies  on  each  wing  was  chosen,  with  the  most 
outboard  body  located  at  the  inboard  edge  of  the  aileron.  Figure  17  shows 
a  sketch  of  the  body  locations  on  the  wing. 


SECTION  III 


WIND  TUNNEL  TEST 

1.  TEST  FACILITY 

The  Arnold  Engineering  Development  Center  Propulsion  Wind  Tunnel  (16T) 
is  a  variable  density,  continuous-flow  closed  circuit  tunnel  capable  of 
being  operated  at  Mach  numbers  from  0.2  to  1.5,  stagnation  pressures  from 
120  to  4000  psfa,  and  Reynolds  numbers  from  0.5  to  5.5  million  per  foot. 

The  test  section  is  16  ft  square  by  40  ft  long  and  is  enclosed  by  60-degree 
inclined-hole  perforated  walls  of  six  percent  porosity. 

2.  MODEL  AND  INSTRUMENTATION 

A  0.044  scale  C-141B  model.  Figure  18,  was  used  in  this  test.  The 
wing  and  empennage  are  constructed  of  4340  steel  and  the  fuselage  is 
aluminum.  The  other  model  components  are  made  from  combinations  of  metal, 
plastic  and  fiberglass.  A  list  of  the  model  components  with  identification 
symbols  is  in  Table  1  and  model  dimensional  data  is  given  in  Table  2.  A 
three  view  sketch  of  the  model  is  shown  in  Figure  19. 

The  forward  12%  of  the  wing,  full  span,  * s  removable.  Modeless  putty 
was  used  to  extend  the  modification  to  17%  chord  near  the  wing  tip.  Two 

leading  edge  configurations  were  investigated  during  this  entry;  the  basic 
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C-141  wing,  W  ,  and  the  2-D  equivalent  approach  wing,  W 

A  total  of  118  static  pressure  orifices  were  located  on  the  right-hand 
wing,  upper  and  lower  surfaces,  at  four  spanwise  stations  as  illustrated  in 
Figure  20  and  listed  in  Table  3.  All  the  orifices  were  connected  to  and 
recorded  from  a  four-module,  48-port  scanivalve  assembly  located  in  the 
model  nose  section.  In  addition,  internal  cavity  pressure  and  differen¬ 
tial  across  the  balance  were  measured  from  two  single  orifices  located  fore 
and  aft  in  the  blade  cavity  and  a  single  orifice  located  aft  of  the  balance. 


The  model  was  mounted  on  a  Task  3500-B  MK-I-3.5  inch  internal  balance 


supported  by  a  blade-sting  system.  This  is  a  six-component,  force  type, 
strain  gage  balance.  Its  design  limits  and  calibration  are  discussed  in 
Reference  10.  A  forward  blade  support  attached  to  the  balance  along 
the  fuselage’s  forward  lower  centerline  was  used  to  minimize  support 
interference  effects  on  the  aft  fuselage.  The  blade  lower  section  was 
attached  to  a  sting  adapter  and  the  AEDC  PWT  16T  "A"  sting. 

The  sting  pitch  and  roll  angles  were  sensed  from  synchrotransmitters . 
The  model  attitude  was  obtained  from  a  combination  of  sting  attitude  and 
sting-balance  deflections  under  aerodynamic  loads.  In  addition,  an  angular 
position  indicator  mounted  in  the  model  was  used  as  a  back-up  pitch  in¬ 
dicator. 

Boundary  layer  transition  was  fixed  on  the  model  by  means  of  Ballotini 
glass  beads  applied  in  strips.  The  transition  strip  location,  width,  and 
bead  diameter  for  each  model  component  are  shown  in  Table  4. 

3.  TEST  CONDITIONS  AND  SCHEDULE 

All  configurations  were  run  at  a  Reynolds  number  of  4.8  million  per 
foot.  (4.69  million  based  on  model  wing  MAC)  and  a  nominal  frees tream 
Mach  number  range  from  M  =  0.60  to  0.83.  Six-component  force  and  pressure 
measurements  were  obtained  at  zero  yaw  over  an  angle-of-attack  range  from 
-4  to  +4  degrees.  Flow  visualization  photos  of  fluorene  sublimation  were 
obtained  to  verify  transition  fixing.  A  summary  of  the  test  program  is 


shown  in  Table  5. 


4. 


DATA  REDUCTION 


Force  and  moment  measurements  were  normalized  to  coefficient  form  in 

the  body  and  stability  axis  system  using  the  following  dimensions: 

2 

Wing  area,  ft  6.247 

Mean  aerodynamic  chord  (MAC),  in  11.724 

Span,  in  84.302 

Reference  Moment  Center,  in  FS  40.605 

WL  10.428 
BL  0.00 

The  balance  output  data  were  reduced  to  coefficient  form  based  on  in¬ 
compressible  dynamic  pressure  and  AEDC's  balance  calibration  prior  to  the 
test.  Blockage  and  tunnel  wall  effects  were  assumed  to  be  negligible  due 
to  the  small  model- tunnel  size  ratio  and  porous  walls  of  the  test  section. 
No  corrections  for  effects  of  blade-sting  tare  and  interference,  nacelle 
internal  drag  or  flow  angularity  were  applied,  as  the  objective  of  this 
test  was  to  identify  drag  increments  between  the  basic  and  modified  con¬ 
figurations  . 

Pressure  data  were  recorded  from  the  +  12.5  psid  transducers  contained 
in  the  48  SGM  scanivalve  module  unit  and  reduced  to  coefficient  form 


where 

C  =  Pressure  coefficient 
P 

P^  -  Model  pressure 

Pg  =  Freestream  or  reference  static  pressure 
q  =  Freestream  dynamic  pressure. 
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During  the  acquisition  of  pressure  data,  computer  evaluation  c 
sure  rate-of -change  was  used  and  the  transducer  output  was  not 
computational  purposes  until  either  the  rate  of  change  was 
acceptable  limits  or  a  maximum  time  delay  was  reached. 
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SECTION  IV 


WIND  TUNNEL  DATA  ANALYSIS 

1.  MODIFIED  LEADING  EDGE  RESULTS 

a.  CHORDWISE  PRESSURE  DISTRIBUTION  COMPARISONS 

Chordwise  pressure  distributions  at  the  four  spanwise  pressure 
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stations  for  the  base,  W  ,  and  modified,  W  ,  wings  are  compared  in  Figures 
21  through  27  for  M  =  0.7,  0.75,  0.77,  0.79,  0.80,  0.81,  and  0.83.  An  angle 
of  attack  of  approximately  1.5°  was  chosen  for  these  comparisons  to  corres¬ 
pond  to  the  C-141  cruise  condition. 

Comparison  of  the  pressure  distributions  at  the  two  inner  stations, 

r)  =  0.193  and  0.418,  at  any  of  the  Mach  numbers  show  very  little  improvement 
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for  the  W  wing  over  the  W  .  There  is  a  slight  increase  in  leading  edge 

suction  but  not  enough  to  give  a  significant  drag  reduction.  At  the  outer 

stations,  n  =  0.637  and  0.793,  there  is  a  definite  increase  in  leading  edge 
37 

suction  for  the  W  wing.  This  increase  in  leading  edge  suction  is  the  type 
of  improvement  sought  and  is  the  major  source  of  drag  reduction  demonstrated 
in  this  test. 

The  leading  edge  modification  had  a  secondary  goal  of  weakening 
the  shock  strength  and  moving  it  aft.  This  never  occurred  experimentally. 
The  only  time  the  shock  was  weakened  was  at  n  8  0.793  and  M  *  0.77;  however, 
the  shock  had  moved  forward  rather  than  aft.  The  shock  moved  aft  only  at 
high  Mach  numbers;  M  =  0.81  and  0.83  for  r|  =  0.793;  and  M  =  0.79,  0.80,  and 
0.81  for  Tl  =  0.637;  however,  the  shock  strength  remained  unchanged.  At 
M  ■  0.75  and  0.77  the  shock  moved  forward.  At  M  =  0.70  a  leading  edge 
shock  formed  on  the  modified  wing  and  this  produced  a  drag  increase. 
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Favorable  effects  on  the  chordwise  pressure  distributions  occurred 
primarily  on  the  outer  half  of  the  wing  with  no  significant  changes  inboard. 
Changes  in  shock  position  were  detrimental  for  M  =  0.77  and  below  and  just 
slightly  favorable  for  M  =  0.79  and  up. 

The  lack  of  shock  movement  on  the  wing  contrasts  sharply  with  the 
shock  movement  found  in  the  2-D  tests.  This  indicates  that  3-D  effects  are 
the  dominating  influence  on  shock  position  and  strength  on  the  wing, 
b.  AERODYNAMIC  CHARACTERISTICS 

(1)  Repeatability 

Repeat  runs  were  made  during  this  test  for  each  configura¬ 
tion.  One  set  of  data  was  taken  near  the  start  of  the  test  of  a  configura¬ 
tion  and  the  second  set  at  the  end  of  the  test.  Figures  28  and  29  show 
these  runs  for  at  M  =  0.60  and  0.77;  Figure  30  shows  at  M  =  0.77. 

The  repeatability  is  excellent  with  the  drag  being  within  +  3  counts 
throughout  the  range.  This  accuracy  of  drag  measurement  was  very  important 
for  this  effort  where  small  increments  between  configurations  are  to  be 
measured.  The  accuracy  of  the  data  acquisition  is  discussed  in  detail  in 
the  report  documenting  the  test  by  Brown  of  AEDC^^. 

(2)  Lift,  Drag,  and  Pitching  Moment  Characteristics 

The  lift,  drag  and  pitching  moment  characteristics  for  all 
four  configurations  tested  are  shown  in  Figures  31  and  32  for  Mach  numbers 
of  0.7  and  0.79.  In  this  section  only  the  and  data  will  be  dis¬ 
cussed.  The  swept  tip  and  anti-drag  body  data  will  be  discussed  in  sections 
IV. 3  and  IV. 4. 

At  M  *  0.7,  Figure  31,  there  are  no  measurable  differences  in 
the  lift  curves  for  the  two  leading  edge  shapes.  The  pitching  moment 
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change,  -0.003,  associated  with  this  change  is  considered  insignificant,  as 

it  represents  less  than  0.1  of  a  degree  of  stabilizer  trim  change  for  the 

C-141.  The  drag  polar,  however,  shows  a  drag  penalty  of  8  counts  at  the 

equivalent  cruise  of  0.5.  This  drag  penalty  is  due  to  the  formation  of 

a  leading  edge  shock  on  the  modified  wing  as  shown  in  Figure  21. 

At  M  =  0.79,  Figure  32,  the  lift  and  pitching  moment  curves 

are  the  same  up  to  the  cruise  C^of  0.50.  However,  the  drag  polar  shows  a 
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20  count  drag  reduction  for  W  at  =  0.50.  This  is  a  drag  reduction  of 
7.4%  at  the  design  cruise  point  and  it  exceeds  the  target  drag  reduction 
of  14  to  16  counts. 

(3)  Drag  Rise  Characteristics 
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Drag  rise  summaries  for  W  and  W  are  shown  in  Figures  33  and 

34  for  a  range  of  lift  coefficients.  These  curves  are  the  result  of  an 

interpolation  of  the  wind  tunnel  data  at  constant  lift  coefficients.  Since 

the  model  was  tested  tail-off  and  untrimmed,  the  cruise  lift  coefficient  is 

approximately  0.5,  corresponding  to  a  trimmed  airplane  lift  coefficient  at 

cruise  of  0,46.  A  comparison  of  the  drag  rise  for  both  wings  at  three  lift 

coefficients  near  the  cruise  value  is  made  in  Figure  35.  A  small  drag 
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penalty  for  W  is  evident  at  speeds  below  M  *  0.75.  At  M  =  0.75  the  drag 

37 

levels  are  identical.  Above  M  =  0.75  the  modified  wing  W  has  a  definite 
drag  advantage.  The  creep  drag  is  reduced  and  drag  divergence  is  delayed, 
c.  INCREMENTAL  DRAG  AND  PITCHING  MOMENT 

Incremental  drag  and  pitching  moment  coefficients  for  the  modified 
wing  are  summarized  in  Figure  36  for  constant  values  of  tail-off  untrimmed 
lift  coefficient.  These  data  are  increments  interpolated  directly  from  the 
test  data.  Scatter,  amounting  to  +  0.0001  in  C^  and  +  0.001  in  C^,  have 
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been  faired  out  in  generating  these  curves.  At  a  typical  cruise  point  for 

the  C-141,  approximately  M  *=  0.77  and  C  =  0.5,  the  W37  modifica- 

LTAIL-0FF 

tion  reduces  drag  by  about  7.5  counts.  This  increases  to  over  20  counts 

of  drag  reduction  at  M  *  0.79.  Pitching  moment  changes  are  less  than  + 

0.004  over  the  cruise  range  and  are  considered  insignificant. 

2.  AERODYNAMIC  CORRELATION  OF  LEADING  EDGE  EFFECTS 

a.  COMPARISON  OF  THEORETICAL  AND  MEASURED  PRESSURE  DISTRIBUTIONS  - 
EQUIVALENT  2-D  METHOD 

In  order  to  compare  the  3-D  test  results  to  the  2-D  theory  the  3-D 
results  must  be  converted  to  2-D  by  means  of  the  relationships  discussed  in 
section  II. 2. b.  Thusv,data  at  the  3-D  design  point  of  M  -  .79  is  compared 
to  the  2-D  theory  at  M  =  .73  and  the  pressure  coefficients  are  ratioed  to 
their  2-D  counterparts  for  comparison.  Figure  37  shows  this  comparison  for 
f|  -  0.418,  0.636  and  0.793. 

At  n  =  0.418,  2-D  viscous  theory  is  shown  and  the  agreement  in 
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shock  position  is  good.  The  leading  edge  pressure  peak  predicted  for  W 
did  not  occur  and  this  may  be  due  to  the  spanwise  influence  of  the  un¬ 
modified  inboard  leading  edge  or  to  the  influence  of  the  fuselage.  The 
discrepancy  in  lower  surfaces  pressures  is  due  to  the  interference  of  the 
wheel  well  fairings  and  pylons  and  nacelles. 

At  f|  *  0.636  and  n  =  0.793  the  agreement  between  theory  and  experiment 
for  both  wings  is  very  good  for  the  leading  edge  pressures.  This  is  the 
critical  area  for  matching  with  the  2-D  equivalent  approach  used  in  the 
design  and  it  confirms  the  soundness  of  the  design  approach.  Experimentally 
the  shock  wave  is  further  forward  than  predicted.  This  is  due  to  the  use 
of  an  inviacid  code  for  the  design. 
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b.  COMPARISON  OF  THEORETICAL  AND  MEASURED  PRESSURE  DISTRIBUTIONS  - 

FLO- 2 2 

The  basic  and  modified  wings  were  analyzed  using  a  3-D  transonic 

wing  code,  FLO-22,  for  purposes  of  comparison.  The  predictions  of  this 

theory  are  compared  to  the  experiment  in  Figure  38.  The  theoretical  pre¬ 
dictions  are  least  accurate  on  the  inboard  wing  since  the  body,  wheel  well 
pods  and  pylons/nacelles  are  not  modeled  in  the  theory.  At  rj  =  0.193  the 
pressure  distributions  do  not  compare  very  well  but  the  shock  position  and 

strength  are  predicted  correctly.  At  n  =  0.418  the  upper  surface  pressures 

match  very  well  for  the  basic  wing,  however,  the  shock  position  is  not  as 
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well  predicted  as  it  is  for  the  rest  of  the  wing.  For  the  W  wing  at 
r|  =  0.418  the  predicted  pressure  peak  does  not  occur  and  the  correlation  of 
upper  surface  pressures  is  not  good.  The  disagreement  in  lower  surface 
Cpfs  for  both  wings  is  due  to  the  interference  of  the  wheel  well  pods  and 
pylon/nacelles.  These  results  underline  the  inadequacy  of  wing  only  codes 
to  design  the  inboard  portion  of  wings.  At  n  =  0.636  and  r)  =  0.793  the 
correlation  between  experiment  and  theory  is  excellent  in  all  aspects.  It 
should  be  noted  that  FLO-22  runs  were  made  at  a  Mach  Number  0.01  above  the 
test  correlations.  Figures  39  and  40  show  this  same  correlation  at 
X]  =  0.793  using  the  0.01  Mach  Number  shift  at  test  Mach  Numbers  of  0.78 
and  0.79. 

Contrasting  Figures  37  and  38  gives  an  indication  of  the  adequacy 
of  using  2-D  inviscid  theory  for  wing  designs.  In  this  case  where  we  were 
only  modifying  the  leading  edge  and  looking  for  changes  in  wave  drag  it 
was  barely  adequate.  However,  more  extensive  modifications  would  call  for 
the  use  of  3-D  codes  which  give  superior  pressure  and  shock  correlations. 
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Figure  38  also  highlights  the  limitation  of  the  wing  alone  code  whenever 
interference  effects  from  the  body,  etc.  are  significant. 

c.  COMPARISON  OF  PREDICTED  AND  MEASURED  DRAG  RESULTS 
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The  incremental  drag  difference  between  W  and  W  experimentally 
is  compared  to  the  theoretical  predictions  in  Figure  41.  This  comparison 
is  at  the  untrimmed  cruise  lift  coefficient  of  0.5.  The  two-  and  three- 
dimensional  theories  agree  well  with  each  other  with  the  most  significant 
difference  being  at  M  =  0.80  which  is  well  into  the  drag  rise.  The  experi¬ 
mental  data  correlates  very  well  for  0.75  £  M  £  0.79  and  this  substantiates 
the  design  technique  used  in  the  effort.  At  M  -  0.80  the  experimental  drag 
increment  is  over  twice  that  predicted  by  theory,  however,  this  is  into  the 
drag  rise  area  where  neither  theory  is  considered  accurate.  At  M  =  0.70 
neither  theory  predicted  the  formation  of  a  leading  edge  shock  wave  which 
produced  the  drag  increase  shown  experimentally. 

3.  SWEPT  TIP  RESULTS 

The  addition  of  the  swept  wing  tip  to  the  modified  wing  has  very  little 
effect  on  the  chordwise  pressure  distributions  as  shown  in  Figure  42  for 
M  =  0.79  and  a  -  1.0  degrees.  At  the  two  outboard  span  stations  there  is 
a  slight  rearward  shift  in  the  shock  which  indicates  a  favorable  effect  on 
the  tip  isobars. 

Figures  31  and  32  compare  the  lift,  pitching  moment  and  drag  character¬ 
istics  tip  on  (triangle  symbol)  to  tip  off  (circle  symbol)  at  M  *  0.70  and 
M  *  0.79.  All  of  these  coefficients  are  based  on  a  wing  area  of  6.247 
square  feet  and  the  added  area  due  to  the  wing  tip  is  ignored.  The  addi¬ 
tion  of  the  wing  tip  increases  the  lift  at  a  given  angle  of  attack  and  as 
Figure  32  indicates  may  delay  the  break  in  the  lift  curve  slope.  The 
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additional  lift  at  a  =  1.0°,  approximately  AC  =  0,015,  is  primarily  con- 

Li 

centrated  at  the  tip  as  the  pressure  data  in  Figure  42  shows.  This  addi¬ 
tional  lift  is  located  behind  the  center  of  gravity  of  the  airplane  and 
creates  the  additional  nose  down  pitching  moment  shown  in  Figures  31  and 
32.  At  =  0.5,  the  pitching  moment  increment  is  AC^  =  -0.025  which 
corresponds  to  approximately  one-half  degree  of  stabilizer  trim  change  for 
the  C-141.  The  expected  large  savings  in  wing  induced  drag  due  to  the  9.5 
percent  increase  in  wing  span  were  realized  as  shown  in  the  drag  polars  in 
Figures  31  and  32. 

The  drag  rise  characteristics  for  the  modified  wing  with  the  swept  tip 
are  shown  in  Figure  43  and  are  compared  in  Figure  44  with  data  for  the 
modified  wing  without  the  swept  tip.  At  C  =0.5  and  for  Mach  numbers 
below  0.76,  the  induced  drag  reduction  amounts  to  about  17  counts.  At 
M  =  0.79  the  drag  reduction  increases  to  25  counts,  indicating  a  reduction 
in  compressibility  drag  of  8  counts.  This  reduction  in  compressibili ty 
drag  indicates  that  the  swept  tips  are  improving  the  shock  losses  over  the 
wing,  especially  at  the  tip  where  the  isobars  become  unswept  on  the  wing 
without  the  swept  tip. 

These  changes  in  aerodynamic  characteristics  are  identical  to  those 
observed  in  the  earlier  test^  where  the  swept  tip  was  added  to  the  basic 
wing  rather  than  the  modified  wing.  This  indicates  that  the  effects  of  the 
leading  edge  and  the  swept  tip  modifications  are  independent  and  additive 
as  we  had  expected  them  to  be. 

Because  the  wing  loading  in  the  region  of  the  tip  is  altered  by  the 
swept  tip  modification,  the  center  of  pressure  for  the  wing  is  shifted  out¬ 
board  a  significant  amount.  In  the  earlier  report^  this  incremental  shift 


in  the  center  of  pressure  at  cruise  was  found  to  be  2.3  percent  of  the  wing 
semi-span.  This  change  in  center  of  pressure  could  be  reduced  or  eliminated 
by  uprigging  the  ailerons. 

4.  ANTI-DRAG  BODY  RESULTS 
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The  anti-drag  bodies,  eight  per  side,  were  added  to  the  W  wing  with 

the  swept  tips  installed  to  evaluate  the  combined  effects  of  all  three 

modifications.  The  earlier  entry had  tested  the  anti-drag  bodies  on  the 
12 

basic  wing,  W  ,  only. 

The  effect  of  the  anti-drag  bodies  on  the  chordwise  pressure  distribu¬ 
tions  is  shown  in  Figure  45.  The  increase  in  aft  loading  and  the  aft  shock 
movement  clearly  indicate  that  the  design  objective  of  increasing  camber 
was  achieved  with  the  anti-drag  bodies.  Station  r|  =  0.793  is  outboard  of 
the  last  anti-drag  body;  however,  the  influence  of  the  anti-drag  bodies  is 
very  evident  on  the  upper  surface  indicating  good  carry-over.  The  lift, 
drag  and  pitching  moment  characteristics  for  this  configuration  are  shown 

by  the  plus  symbol  on  Figures  31  and  32.  This  data  can  be  compared  to  the 
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W  plus  swept  tip  data,  triangle  symbols,  to  get  the  incremental  effect 
due  to  the  addition  of  the  anti-drag  bodies.  The  angle  of  attack  needed 
to  produce  a  cruise  lift  coefficient  is  almost  a  degree  less  due  to  the 
camber  effect.  The  nose  down  pitching  moment  increases  by  -0.015  at 

=  0.5  and  this  would  require  approximately  0.3  degree  of  stabilizer  trim 
change.  The  drag  polar s  show  a  definite  drag  penalty  at  low  lift  co- 
efficients  and  a  slight  drag  reduction  at  lift  coefficients  above  0.5. 

The  drag  rise  characteristics  for  the  anti-drag  bodies  are  shown  in 
Figure  46.  Figure  47  shows  the  effect  of  the  anti-drag  bodies  on  the  drag 
rise  of  the  modified  wing  with  swept  tip.  This  modification  only  provides 
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a  drag  reduction  above  M  =  0.79  at  the  design  cruise  lift  of  C L  *  0.5. 

These  drag  results  are  not  quite  as  good  as  the  results ^  when  the 
anti-drag  bodies  were  tested  on  the  basic  wing.  Thus  the  drag  effects  of 
the  anti-drag  bodies  do  not  appear  to  be  totally  additive  to  the  effects 
of  the  leading  edge  and  swept  tip  modifications. 


SECTION  V 


PERFORMANCE 

1.  SUMMARY  OF  EFFECTS  ON  C-141B  CRUISE  PERFORMANCE,  PRODUCTIVITY  AND 
MISSION  FUEL  FOR  THE  MODIFIED  LEADING  EDGE. 

The  baseline  configuration  for  the  performance  evaluation  in  this 
study  is  the  C-141B  aircraft  as  flight  tested  in  FY80.  Performance  im¬ 
provements  for  the  leading  edge  modification  were  added  incrementally  to 

the  performance  for  this  aircraft.  No  corrections  were  applied  to  the 

12  37 

tested  drag  increments  between  the  W  and  W  configurations. 

The  cruise  lift  coefficient  for  the  C-141  is  typically  0.46,  trimmed, 
and  this  has  been  chosen  for  the  cruise  summary  comparisons.  This  value 
corresponds  to  approximately  C  =0.5  for  the  tail-off  un trimmed  test 

Lj 

results  reported  herein.  Figure  48  summarizes  the  effect  of  the  leading 
edge  modification  on  the  C-141B  drag  rise  and  range  parameter,  M(L/D) ,  for 
this  cruise  condition.  Two  conclusions  can  be  drawn  from  these  results. 
First,  an  increase  in  the  range  parameter  is  indicated  above  M  =  0.74. 
Second,  the  drag  decrease  Improves  substantially  as  the  Mach  number  is  in¬ 
creased  to  the  design  point,  M  =  0.79. 

The  performance  curve  for  the  leading  edge  modification  only  was  put 
into  the  C-141B  performance  deck  to  determine  relative  specific  range  values 
for  the  basic  and  modified  configurations.  The  differences  in  cruise  speeds 
resulted  in  changes  in  cruise  altitude  and  specific  fuel  consumption  (SFC) 
which  offset  the  fuel  savings  due  to  the  drag  reductions  obtained  so  that 
net  mission  fuel  remained  the  same.  This  result  was  very  disappointing 
since  from  a  fuel  savings  standpoint,  the  increased  cruise  speed  was  not 
beneficial.  However,  the  increased  cruise  speed  could  be  utilized  to 


increase  the  productivity  of  the  C-141  fleet. 

2.  SUMMARY  OF  EFFECTS  OF  THE  SWEPT  WING  TIPS  ON  C-141B  CRUISE  PERFORMANCE. 

The  effect  of  the  swept  wing  tips  added  to  the  C-141B  with  a  modified 
leading  edge  is  also  shown  in  Figure  48.  The  effects  of  these  two  modifica¬ 
tions  have  been  shown  to  be  additive,  thus  the  effect  of  the  swept  tips  can 
be  extracted  from  the  M(L/D)  vs.  M  data  shown.  The  swept  tip  offers  a  sub¬ 
stantial  improvement  in  M(L/D)  over  the  entire  Mach  number  range  tested, 

0. 70  <  M  0 . 83  and  increases  the  Mach  number  for  best  M(L/D)  to  M  =  0.79. 
This  can  be  attributed  primarily  to  the  9.5  percent  increase  in  wing  span. 
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SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 

An  analytical  and  wind  tunnel  test  program  has  been  completed  on  a 
leading  edge  modification  to  the  C-141B  wing  in  order  to  evaluate  cruise 
drag  reduction.  Two  additional  drag  reduction  concepts,  swept  wing  tips 
and  anti-drag  bodies,  were  tested.  The  principal  conclusions  follow. 

a.  A  successful  application  of  a  low  cost  2-D  transonic  numerical 
optimization  method  (CONMIN)  to  design  an  improved  leading  edge  contour  for 
the  C-141B  wing  has  been  demonstrated.  For  this  unique  problem  where  only 
the  forward  17  percent  of  the  wing  was  modified  the  use  of  a  low  cost  2-D 
code  instead  of  the  more  costly  3-D  code  was  found  to  be  adequate  for  pre¬ 
diction  of  leading  edge  pressures  but  not  to  be  adequate  for  shock  pre¬ 
diction.  The  correlation  of  leading  edge  pressures  on  the  outboard  half 

of  the  wing  was  excellent.  The  correlation  at  r\  =  0.418  was  not  good; 
however,  this  station  is  between  the  engines  and  even  current  3-D  codes 
would  not  be  expected  to  give  accurate  results  here. 

b.  The  advantages  of  a  3-D  transonic  wing  code  over  the  2-D  code  is 
shown  in  the  accurate  prediction  of  shock  strength  and  location  at  all  four 
span  stations.  The  correlation  of  Cpfs  are  excellent  for  the  outboard  span 
stations . 

c.  The  need  for  body  and  pylon/nacelle  modeling  is  shown  by  the  lack 
of  pressure  correlation  at  the  two  inboard  span  stations.  The  differences 
between  the  wing  code  predictions  and  experiment  can  be  attributed  directly 
to  aerodynamic  interference  from  these  components. 

d.  The  creation  of  a  leading  edge  pressure  peak  did  reduce  the  creep 
drag  for  this  aircraft.  Comparison  of  the  experimental  pressure  distribu¬ 
tions  reveals  that  the  change  in  leading  edge  pressures  is  the  only 
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significant  difference  between  the  two  wings;  therefore,  this  must  be  the 
reason  for  the  drag  reduction* 

e.  The  increased  cruise  speed  capability  provides  a  potential  benefit 
with  respect  to  productivity  of  the  C-141.  By  taking  advantage  of  the 
higher  cruise  speed,  mission  time  is  decreased,  which  either  reduces  flight 
hours  for  the  fleet  or  increases  the  total  fleet  productivity  capability. 

f.  Substantial  drag  reduction  was  measured  at  the  design  Mach  number 
of  0.79.  This  resulted  in  increases  in  aerodynamic  cruise  efficiency, 
M(L/D) ,  and  in  the  Mach  number  for  maximum  M(L/D) .  However  increases  in 
engine  SFC  as  Mach  number  increases  negate  the  expected  fuel  savings. 

g.  The  effect  of  the  swept  wing  tip  was  confirmed  as  being  additive 
to  the  effects  of  the  leading  edge  modification  since  the  changes  due  to 
adding  the  swept  tip  to  the  modified  wing  were  virtually  identical  to  the 
changes  when  it  was  added  to  the  basic  wing. 

h.  The  effect  of  adding  the  anti-drag  bodies  to  the  modified  wing 
with  swept  tips  was  not  as  great  as  the  effect  of  adding  them  to  Che 
basic  wing.  This  was  expected  since  the  design  of  the  anti-drag  bodies  was 
based  on  the  basic  airfoil  and  wing  shapes. 

2.  In  view  of  the  above  discussion  and  conclusions,  the  following  is 
recommended. 

a.  Future  work  on  leading  edge  modifications  should  address  the 
cruise  speed  range  of  0.74  <_  M  j<  0.76  to  avoid  detrimental  changes  in 
engine  SFC. 

b.  A  different  design  objective  should  be  chosen  for  the  CONMIN 
optimizer.  This  effort  used  minimum  wave  drag  as  the  design  objective. 
Other  quantities  such  as  total  drag  or  leading  edge  suction  could  be 
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utilized. 


c.  Active  consideration  of  low  speed  performance  (take-off,  landing, 
and  air-drop)  should  be  part  of  any  future  leading  edge  design  since  this 
wing  is  the  result  of  considerable  compromise  between  low  speed  and  cruise 
requirements. 

d.  A  3-D  wing/body  code  should  be  used  for  redesign  of  the  in-board 
half  of  the  wing.  At  40%  semi-span  the  need  for  pylon/nacelle  modeling  is 
very  evident. 

e.  The  swept  wing  tip  should  be  studied  further.  A  load  alleviation 
scheme  using  up  rigged  ailerons  should  be  demonstrated.  These  results 
should  then  be  compared  to  a  winglet  designed  for  the  C-141B  and  the  rela¬ 
tive  benefits  and  penalties  should  be  clearly  identified. 


TABLE  1 

CONFIGURATION  SYMBOLS 


SYMBOLS 


,12 


K 


19 


N 


8 


W 


12 


W 


37 


Z 


fl 


Z 


G21 


COMPONENT 

FUSELAGE  -  C-141B  Base;  Basic  C-141A 
Fuselage  with  280  inch  Full  Scale 
Extensions 

PYLONS  -  Engine  Nacelles 

NACELLES  -  Flow  Through,  with  Inlet 
Spinners 

WING  -  Basic  High/Low  Speed  W^ 

Steel  Wing  with  Removable  Leading 
Edges  (Forward  12%  C^) .  3  Panels  L.E. 
Each  Semi-Span 

12C 

WING  -  W  With  2-D  Equivalent 
Approach  Leading  Edge 

ANTI-DRAG  BODIES  -  Flap  Track  Fairing 
Type,  Wing  Trailing  Edge;  8  per 
Semi-Span 

WHEEL  WELL  FAIRING 


Z 


t6 


WING  TIP  -  Swept;  Chord/Span  Ratio  - 
1.33/ i 


ZW7A 


WING-FUSELAGE  FILLET  -  C141A  Pro- 
duction  Fillets 


TABLE  2 

MODEL  DIMENSIONAL  DATA 


FUSELAGE  -  B12 

Length,  Inches 

82.350 

2 

Max.  Frontal  Area,  Ft 

0.305 

Max.  Equiv.  Diameter,  Inches 

7.480 

Fineness  Ratio  (L/^) 

11.009 

Nose  Location,  FS 

3.098 

Fuselage  Reference  Line  (FRL) ,  WL 

8.800 

PYLONS  -  K19 

INBOARD 

OUTBOARD. 

2 

Area,  Ft 

0.090 

0.093 

Span ,  Inches 

1.453 

1.517 

MAC ,  Inches 

8.800 

8.800 

Thickness,  Streamwise,  %  Chord 

8.000 

8.000 

Sweep  of  Leading  Edge,  Degrees 

73.000 

73.000 

NACELLES  -  N8 

2 

Area,  Side  per  Nacelle,  Ft 

0.168 

Length,  Inches 

8.256 

Max.  External  Diameter,  Inches 

2.900 

Internal  Diameter,  Inches  -  Inlet 

2.244 

-  Exit 

1.940 

2 

Internal  Area,  In  -  Inlet 

3.954 

-  Exit 

2.956 

Fineness  Ratio,  (^/d)  Ext. 

2.843 

Toe-In,  Degrees  -  Inboard 

2.0 

-  Outboard 

1.0 

Location  of  Inlet  -  INBOARD,  FS 

29.104 

WL 

8.479 

BL 

12.272 

-  OUTBOARD  FS 

32.709 

WL 

8.163 

BL 

20.103 
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TABLE  2  (CONT ’D) 


WING  -  W12 

Area,  Ft2  (S)  6.247 

Span  (Equiv.),  Inches/Ft(b)  84.302/7.025 

Aspect  Ratio  (A)  7.9 

Taper  Ratio  (  X  )  0.373 

Thickness  Ratio  (t/c)  -  Root  (BL  0.00)  0.130 

-  Inboard  Break  (BL  17.804)  0.112 

-  Outboard  Break  (BL  18.778)  0.110 

-  Tip  0.100 

-  Mean  0.113 

Sweep  of  25%  chord.  Degrees  -  Inboard  23.734 

-  Outboard  25.025 

Sweep  of  Leading  Edge,  Degrees 

-  Inboard  28.253 

-  Outboard  27.285 

Dihedral,  Degrees  -  Inboard  -0.941 

-  Outboard  -1.195 

Incidence,  Degrees  from  FRL  @  BL  0.00  (i)  4.891 

Twist,  Degrees  -  Root  (BL  0.00)  0.00 

Inboard  Break  (BL  17.804)  -2.201 

Outboard  Break  (BL  18.778)  -2.279 

Tip  (BL  42.151)  -5.584 

Chord  Lengths  (Projected),  in.-  Root  (c^)  BL  0.0  17.547 

MAC  (c)  BL  17.309  11.724 

Inbd.  Brk.  BL  17.804  10.591 

Outbd.  Brk.  BL  18.778  10.379 

Tip  (ct)  BL  42.151  5.803 

Chord  Locations  -  Root  (L.E.)  FS  28.545 

WL  13.045 


BL 

0.00 

(MAC  (25%) 

FS 

40.605 

WL 

11.629 

BL 

17.309 

TIP(L.E.) 

FS 

50.531 

WL 

10.219 

BL 

42.151 

*NOTE :  Wind  Tunnel  data  based  on  projected  wing  planform  with  root 
chord  @  0°  incidence  and  aero  L.E.  and  100%  chord  T.E. 

Aero  data  analysis  based  on  area  measured  from  L.E.  to  T.E. , 
and  25%  MAC  location  referenced  to  wing  reference  plane. 

W.T.  Data  25%  MAC  location  =  24.1%  MAC  Aero  Analysis  Data. 

p  on 

WHEEL  WELL  FAIRINGS  -  Z 


Length,  Inches 

17.776 

Max.  Frontal  Area  Per 

Side,  Inches 

7.480 

Max.  Equiv.  Diameter, 

Inches 

3.086 

Fineness  Ratio 

5.760 

Leading  Edge  Location, 

FS 

34.936 
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TABLE  3 

WING  PRESSURE  ORIFICE  LOCATIONS 


I 


BL  8. 

127 

BL  17.634 

BL  26. 

825 

BL  33. 

437 

n=  C. 

193 

*1  = 

0.418 

V  =  0. 

637 

n=  0. 

793 

(X/C's  FOR  ORIFICE  LOCATIONS ) 

UPPER 

DOWER 

UPPER 

LCWER 

UPPER 

LOWER 

UPPER 

LOWER 

0.0 

0.02 

0.0 

0.02 

0.0 

0.02 

0.0 

0.02 

0.015 

0.05 

0.015 

0.05 

0.015 

0.05 

0.015 

0.05 

0.03 

0.10 

0.03 

0.10 

0.03 

0.10 

0.03 

0.10 

0.05 

0.20 

0.05 

0.15 

0.05 

0.15 

0.05 

0.15 

0.07 

0.30 

0.07 

0.20 

0.07 

0.20 

0.07 

0.20 

0.09 

0.40 

0.09 

0.30 

0.09 

0.30 

0.09 

0.30 

0.11 

0.50 

0.11 

0.40 

0.11 

0.40 

0.11 

0.40 

0.20 

0.63 

0.15 

0.50 

0.15 

0.50 

0.15 

0.50 

0.30 

0.80 

0.20 

0.65 

0.20 

0.65 

0.20 

0.65 

0.40 

0.95 

0.25 

0.85 

0.25 

0.85 

0.25 

0.85 

0.50 

0.30 

0.30 

0.30 

0.63 

0.37 

0.35 

0.35 

0.80 

0.40 

0.40 

0.40 

0.95 

0.45 

0.45 

0.45 

1.0 

0.50 

0.50 

0.50 

0.60 

0.60 

0.60 

0.70 

0.70 

0.70 

0.80 

0.80 

0.80 

0.90 

0.90 

0.90 

0.95 

0.95 

0.95 

1.00 

1.00 

1.00 

37 


Table  4 


Boundary- Layer  Transition  Strip  Location 


Component 

Location 
(In.  from 
Leading  Edge) 

width 

(In.) 

Bead  Diam. 

(In. ) 

Wing  Upper  Surface 

0.  70 

mmm 

0.0031 

Wing  Lower  Surface 

0.70 

1 

0.0031 

Fuselage 

0.90 

0.  10 

0.0045 

Pylons 

0.20 

0.  05 

0.0038 

Nacelles  (External) 

O 

00 

« 

o 

0.05 

0. 0038 

Notes:  1)  All  measurements  along  body  x-axis. 

2)  After  RN  315  the  bead  diameter  of  the  wing 
upper  surface  was  increased  to  0.0031. 
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TABLE  5 

TEST  PROGRAM  SUMMARY 


a)  Model  Configurations 


RN  No. 
START 

RN  No. 
END 

CONF. 

CONFIGURATION 

315 

342 

1 

B12, 

W37, 

ZW7,  ZG21,  K19,  N8, 

ZT6 ,  ZF1 

347 

352 

3 

B12, 

W37 , 

ZW7,  ZG21 ,  K19,  N8, 

ZT6 

359 

368 

5 

B12 , 

W37 , 

ZW7,  ZG21 ,  K19 ,  N8 

381 

398 

10 

B12, 

VJ12 , 

ZVJ7 ,  ZG21 ,  K19 ,  N8 

b)  Test  Program 


Description  & 

Configuration 

PWT  RN 

@  Indicated  Mach  Number 

.600 

.700 

.750 

.770 

.780 

.79 

.80 

.81 

.830 

I .  SUBLIMATION 

■ 

Conf  1 

I 

315 

■ 

317 

II.  Force  &  Pressure 

■ 

Conf  1 

338 

339 

340 

341 

342 

||jl 

322 

Conf  3 

347 

348 

349 

350 

351 

352 

Conf  5 

359 

360 

361 

368 

363 

364 

365 

366 

367 

362 

1 

Conf  10 

381 

384 

i 

385 

1 

392 

387 

388 

389 

390 

391 

383 

i 

386 

i 

| 

i 

III.  SUBLIMATION 

I 

Conf  10 

i 

i 

1 

395 

i 

1 

I 

1 

398 
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IMPORTANCE  OF  LEADING  EDGE  REGION 


MACH  NUMBER 

Research  Airfoils-Example  of  the  Effect  of  Leading  Edge  Region  on  Airfoil  Drag  Rise 


SYM  LE 
©  BASIC 

0  3A 

<S>  4 

A  6 


X/C 


b.  Comparison  of  Experimental  Pressure  Distributions 


Figure  3.  Continued 
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ESTIMATED  C-141BWITH  "A"  FILLET 


DESIGN  ORDINATES 
MEASURED  MODEL  ORDINATES 


Fiaure  5.  Comparison  of  W  and  W  Measured  and  Design  Ordinates  at  *0.793 


Effect  of  Leading  Edge  Change  on  Airfoil 


COUNTS 


BASIC  (W 


Figure  14.  Comparison  of  Sparwise  Variation  of  Wave  Drag 
Increments  for  the  Modified  Leading  Edges 
Relative  to  the  Baseline 


LOWER  SURFACE  MODIFICATION 


M  =  0.74  Cj  *  0.60  7) “  0 .389 

- BASIC  AIRFOIL 

a  =  2. 7°,  C^  =.0024 


W 


b.  Effect  on  Theoretical  Pressure  Distribution 
Figure  16  .  C~I4I  Airfoil  Aft  Camber  Change 


***** 


Figure  IB.  Photograph  of  the  0.044  Scale  C-141B  Model 
Installed  at  AEPC  1 Bt 


FS=40 .605 


Figure  19.  Three  View  Sketch  of  the  0.044  Scale  C-141B  Model 


FS-40.605 


SYM 

CONF 

PART 

ALPHA 

□ 

w12 

386 

1.59 

O 

w37 

362 

1.56 

b .  n  «  0.418 
Figure  23.  Continued 
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SYM  CONF  PART  ALPHA 

□  W12  386  1.59 

O  W37  362  1.56 


d.  n  “  0.793 
Figure  23.  Concluded 
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Figure  24.  Effect  of  W  Leading  Edge  Modification 
on  Chordwise  Pressure  Distributions  at 
M  *  0.79 


MICROCOPY  NlMH.ll  l  ION  H  SI  CM  WO 


SYM 


CONF 


PART 


□ 

w12 

390 

o 

w37 

366 

T 


T 


~T 

n  r 


n 

o. 


a.  n  -  0.193 


Figure  27 .  Effect  of  W*^  Leading  Edge  Modif 
Chordwise  Pressure  Distributions 
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CONF 


PART 


ALPHA 


a.  Lift  and  Pitching  Moment 
Pitching  Moment  Characteristics 


Figure  31.  Concluded 


MACH  NUMBER 


b.  Drag  Increments 


Figure  36.  Concluded 


□  □ 


□ _ □ 


'I'  0.0 


o  J  0A  0-* 

°f  X/C 


a.  n  =  0.418 

Figure  37.  Comparison  of  2-D  Equivalent  Theory  to 
Experimental  Pressure  Distributions  at 

*rest  -  °-79 


110 


o 


Figure  39,  Comparison  of  FLO-22  Predictions  to  Experimental 
Pressure  Distributions 
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SYM  MACH 

FLO  22  -  0.80 

TEST  □  0.79 


ALPHA 

1.00 

1.10 


rit  Ji 

LJ 


-  FLO  22 

-  BAUER  IN VISCID 

O  TEST  RESULTS 
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Figure  41.  Comparison  of  Predicted  and 

Measured  Drag  Increments  for  Vr 
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SYM 

CONF 

PART 

ALPHA 

□ 

w37 

364 

1.06 

o 

37 

W  +TIP 

350 

0.94 

b.  n  -  0.418 


42.  Continued 
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SYM  CONF  PART  ALPHA 

□  W37  364  1.06 

A  W37+TIP  350  0.94 


d.  n  =  0.793 
re  42.  Concluded 


123 


SYM 


CONF 


PART 


ALPHA 


□ 

w37 

350 

0.94 

o 

37 

W  4-TIP 
4ADB 

340 

1.01 

a.  n  ■  0-193 

Effect  of  Anti-Drag  Bodies  on  Chordwise  Pressure 
Distributions  at  M  *  0-79 


126 


SYM 

CONF 

PART 

ALPHA 

□ 

37 

W  '+TIP 

350 

0.94 

A 

37 

W  +TIP 
+ADB 

340 

1.01 

d.  o  -  0.793 


45.  Concluded 
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